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Abstract 

We study the Andreev reflection (AR) in a three terminal mesoscopic hybrid 
system, in which two ferromagnets (Fi and F2) are coupled to a supercon- 
ductor (S) through a quantum dot (QD). By using non-equilibrium Green 
function, we derive a general current formula which allows arbitrary spin po- 
larizations, magnetization orientations and bias voltages in Fi and F2. The 
formula is applied to study both zero bias conductance and finite bias current. 
The current conducted by crossed AR involving Fi, F2 and S is particularly 
unusual, in which an electron with spin a incident from one of the ferromag- 
nets picks up another electron with spin a from the other one, both enter S 
and form a Cooper pair. Several special cases are investigated to reveal the 
properties of AR in this system. 
PACS numbers: 74.50.+r, 73.40.Gk, 75.70.Pa 
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I. INTRODUCTION 



Electron has spin as well as charge. The application of the electron spin property opens 
a fruitful field in the transport of ferromagnetic materials, such as the discovery of GMR 
and TMR effects On the other hand, there are growing interests on the mesoscopic 
normal-metal / superconductor (N/S) hybrid system 0, in which Andreev reflection (AR) 
at N/S interface plays an important role in low bias voltage regime [0]. In AR process, an 
electron incident with energy E and spin a picks up another electron with energy -E and 
spin a, both enter S and form a Cooper pair, leaving a Andreev reflected hole in N side. One 
may expect that the interplay of the spin property of AR process and the spin-dependent 
transport in ferromagnetic materials will add new physics to mesoscopic hybrid systems, 
and to the future applications of spintronics. 

Several works have been devoted to this issue. In the pioneering work of de Jong et al. , 
the transport of a ferromagnet / superconductor (F/S) junction was studied by scattering 
matrix formalism. The conductance of AR is shown to be strongly affected by the spin 
polarization of F. The idea was verified by recent experiments in F/S thin film nanocontact 
and F/S metallic point contact p. Especially, in Ref. 0, Soulen et al. successfully 
determined the spin polarization at the Fermi energy for several metals by measuring the 
differential conductance of F/S metallic point contact. Further calculations [0 implied that 
the Fermi velocity mismatch between F and S also affects AR conductance of F/S contact, 
and the conductance may even be enhanced in presence of spin polarization. In addition to 
simple F/S junction, F/S contact with S in d-wave symmetry 0,^, F/S nanostructure with 
giant proximity effect P,pi]|, and more complicated structures such as FSF double junctions 
TTHl3|, SFS double junctions [1T4|-|T^, S/F superlattices 0], (NF)„S multilayer structures 



T8| , p!9[| were also investigated. 

In this paper, we propose an idea that two sources of spin polarized electrons with 
different orientations are injected into a superconductor, which can be achieved by a three 
terminal mesoscopic F/S hybrid structure shown in Fig.l. In this system, a central quantum 
dot (QD) is coupled via tunnel barriers to two ferromagnetic electrodes (Fi and F2) and a 
superconducting electrode (S) ( hereafter, the system is simply referred as to (Fi,F2)-QD-S ). 
Fi and F2 are assumed to have arbitrary magnetization orientations, spin polarizations, and 
bias voltages. The bias voltage of S is set to zero as the ground. QD is designed to provide a 
link between Fi, F2 and S, so that AR can take place through discrete energy states of QD. 
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Consider the special case that Fi and F2 are fully polarized, AR only involving Fi and S or 
only involving F2 and S are completely suppressed, while the crossed AR involving Fi, F2, 
and S depends strongly on the magnetization orientations of Fi and F2, being suppressed 
if they are in ferromagnetic alignment, enhanced in anti- ferromagnetic alignment. In this 
paper, we will derive a current formula by using non-equilibrium Green function method, 
and investigate several special cases to illustrate the properties of ARs in this system. 

During the preparation of this paper, we are aware that in the recent publication of 
Deutscher et al. |2^, a device consisting of two point contacts between two ferromagnetic 
tips and a superconductor was proposed. For the two tips with fully but opposite spin 
polarizations, they suggested that "mixed" Cooper pair made of electrons coming one from 
each tip can be injected into the superconductor, leading to unusual properties of such 
device. Sec. IV of this paper is partially stimulated by their work. 

The rest of this paper is organized as follows: In Sec. II, we present the model Hamiltonian 
and derive a general current formula for the hybrid system (Fi F2)-QD-S, by non-equilibrium 
Green function method. In Sec. Ill, we study the zero bias conductance, assuming Vi = V2 = 
O"*". The explicit forms of the conductance are presented and numerically studied. In Sec. IV, 
we study the finite bias current with Fi and F2 in anti-ferromagnetic alignment, and the 
fully spin polarized case is discussed in detail. Finally, a brief summary is given in Sec.V. 

II. MODEL AND FORMULATION 

The system under consideration can be described by the following Hamiltonian: 

H = H, + H2 + Haot + H, + HT , (1) 
Hi = ^(cfc - ahi - ^ii)a\^aka , 

ka 

ka' 

Hdot = Eq"^ c\c„ , 

Hs = ^ ^pdljpa + [^^T^-Pi + ' 

pa p 

Ht = Y ^laal^Ca + /i.cj + ^ |^t2f7&Lca + /i-cj + ^ [tsdl^Ca + h.c] . 

ka ka pa 
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Hi and H2 are the Hamiltonians of Fi and F2 in the mean field approximation, with different 
magnetization orientations and chemical potentials. The spin bands of Fi (F2) are split by 
2hi (2/12) due to the exchange energy. The magnetization orientation of Fi is set as z axis, 
while the orientation of F2 as z' axis which has an angle 9 with respect to z axis. The 
operators with the spin- quantization axis z and the operators with the spin-quantization 
axis z' are related by Ds matrix as 



( 



\ 



cos I — sm I 



sm 



cos ■ 



(2) 
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a dot describes the quantum dot, in which only one spin degenerate level is considered and 
the intra-dot interaction is ignored for simplicity. Hg is the Hamiltonian for a BCS super- 
conductor with the chemical potential fixed to zero as the ground. Ht depicts the tunneling 
between QD and Fi, F2 and S, coupling different parts of the system together. 

Since the current through QD can be expressed in terms of the Green functions of QD, we 
first derive the retarded and distribution Green functions by Dyson equation and Keldysh 
equation. To include the physics of Andreev refiections and the spin flip processes in a 
unified formulation, we introduce a 4x4 matrix representation, in which the Green function 
is defined as 



ct cJ ci c| 



)) ■ 



(3) 



Let denote the Fourier transformed retarded Green function of QD , and C can be 
solved by Dyson equation: 

G"= g'+g'-S^C , (4) 

in which g'' is the retarded Green function of an isolated QD and S'' is the self-energy due 
to couplings between QD and leads, g^ can be easily obtained as: 



g' 



1 










1 



a;-|-£;o+iO+ 









1 







-£0+10+ 




^ 






(5) 
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while S'' consists of three parts, 



(6) 



is the self-energy from the coupling between QD and Fi, given by 



\ 



Tit 














Tu 














Tu 

















(7) 



in which and Fij^ are the spin-dependent coupling strengths defined by Fig. = 2t^Ni^ \ti^\ , 
with Ni„ being the density of states of spin a band of Fi. is the self-energy from the 
coupling between QD and F2, given by 

\ 



^ c^Fst + s'^V^i sc(F2t - V21) 



^2 = 







cT2i + s^Fat 











sc(F2t - F2i) 



sc(F2t-F2j c2F2i + s2F2t 



V 







sc(F2T - F2J 







c2F2t + sT 



(8) 



21 J 



in which s = sin | , c = cos | , F21 and F2j^ are defined similarly to Fi-f and Fij^. is the 
self-energy from by the coupling between QD and S, given by 



/ 1 - 



--F,p(a;) 







-^100 



1^ 



V 



0^1 



(9) 



in which Fg = 2nNs \ ts\ , with A^^ being the density of states when the superconductor is in 
normal state, p{uj) is the modified BCS density of states defined by 



1^1 



i^/A' 



\uj\ > A 
\uj\ < A 



(10) 



Thus, C can be obtained by solving Dyson equation, Eq.(4). 
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Let denote the Fourier transformed distribution Green function of QD, and G^ can 
be obtained by Keldysh equation: 



G<= G"S<G" 



(11) 



Notice that the advanced Green function and self-energy are the Hermitian conjugations of 
the corresponding retarded Green function and self-energy. And can be obtained by 
applying the fluctuation-dissipation theorem to each of , S2 and Hf, 



in which 



S< = 



Fi 



E< + S< + E< 
Fi(S?-SO , 

F2(S5-SD , 
F.(S:-SD , 



^ fi ^ 

A 

/i 

A y 



(12) 



(13) 



Fo = 



^ f2 ^ 

A 

/2 

A y 



(14) 



F. = 



O ^ 



0/00 
0/0 

/y 



(15) 



where /i, /i, /s, A and / denote /(a; - \/i), /(a; Vi), f{uj - V2), f{uj + V2) and /(a;), 
respectively, in which /(a;) is the Fermi distribution function. 
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Then, the current flowing from Fi to the QD can be expressed in terms of and G* 



as 



h = /it + In = ^jdu; [(GSi)< + h.c] 



11+33 



(16) 

+ 



in which we have used the compact notations [AB] = A'B'^+A^B" and [ = [ m 

[ Jgg. After some algebra manipulations, the current can be divided into contributions from 
four conducting processes: 

h = lj duj [An(/i - /i) + Ai2(/i - h) + gi.(/i - /.) + gi2(/i - /2)] , (17) 



in which 



^11 - riT(ri-f 1^12!^ + |G32|^) + ri-f(ri-f |Gi4|^ + |G34|^) 



represents the Andreev reflection through Fi-QD-S, 

^12 = {c^r^i + s'r2T)(FiT \G\^\' + Fi^ |G^2l') + 

(c^F2T + s'F2i)(F^|G^4l' + ru|G^4l') + 
5c(F2T - F2T)2Re(F^Gl2Gi: + V.^Gl^GH) 

represents the crossed Andreev reflection through (Fi,F2)-QD-S, 

A A 
2Re( G'G'* + -( 



(18) 



(19) 



Qis — Ti-fF^p 
TuFsp 



\G\, 




1^12 




1^13 






1^31 




1^32 




1*^33 




1 z^** 1 
l*^34l 



+ (20) 



A 
— ( 



"^31"^ 32 



A 
— ( 

UJ 



"^33*-^ 34 J 



represents the single particle tunneling through Fi-QD-S, and p{uj) = ~ ^) is 

the ordinary BCS density of states, 

Q12 = (c^r2T + s'r2i)(F^ \G\,f + Fi^ \Gl,\^) + (21) 

{c'T,^ + s'T,,){T,,\G\,\' + Tn\Gl,\') + 
sc(F2T - F2J2 Re(F^GliG- + T.^Gl^G^*^) 

represents the single-particle tunneling through F1-QD-F2. Similarly, one can obtain the 
current flowing from F2 into QD simply by exchange the index 1 and 2. 

The current formula Eq.(17) is the central result of this work, which can be applied to 
ferromagnetic electrodes Fi and F2 with arbitrary spin polarizations, magnetization orien- 
tations and bias voltages. 



In the following numerical studies, we assume that \eVi\ , \eV2\ < A and ksT <^ A. Qis 
process will vanish because of the factor p and the Fermi function difference (/i — fs). Qu 
process will be ruled out in two special cases: Fi and F2 are either equally biased (Section 
III) or fully but oppositely polarized (Section IV). We will concentrate on AR processes An 
(direct AR through Fi-QD-S) and Au (crossed AR through (Fi,F2)-QD-S ), and investigate 
several special cases to illustrate the properties of these two ARs. 



III. ZERO BIAS CONDUCTANCE 



In this section, we study the zero bias conductance by taking Vi — V2 — 0'^. Since no bias 
voltage between Fi and F2, there is no net single-particle current flowing between them. For 
ksT <^ A, the single-particle current from Fi or F2 to S is also neghgible. Therefore, only 
ARs contribute to the conductance. For simphcity, we set ksT — and £^0 — 0, introduce the 
spin polarization Pp = p^|_^_p|^| , and the spin-averaged coupling strength = ^(F^^ + Tp^), 
with f3 = 1,2 for Fiand F2 respectively. 

First consider the simplest case in which r2 = 0, Fi = F^, Pi = P, Tg = Tn, then 
the three-terminal system (Fi,F2)-QD-S reduces to a two-terminal system F-QD-S, and the 
conductance is easily obtained from the current formula as 

^^^^ - ~h{l-P^ + rr ' ^ ^ 

where r = T^/Tl is the ratio of the two coupling strengths. Analogous to the matching 
condition of the Fermi velocities in F/S contact (i.e., kpikpi = kg), here + = 1 (i.e., 
^L^^Li = r|;) plays the similar role. For r > 1, the matching condition can never be 
satisfied, so Gpos decreases monotonously with the increase of P (Fig. 2a). While for r < 1 
, it exists a certain value of P, say Pq, satisfying Pq + r'^ — 1, so Gpos first increases with 
P , reaches its maximum value ^ at P = Pq, then drops to when P approaches to 1 
(Fig. 2b). This result warns us to be careful to deduce the spin polarization of F from AR 
conductance of F-QD-S. 

Next, consider the general case of the three terminal system (Fi,F2)-QD-S. Similar to 
the composition of polarized light, the total current (or total conductance) of Fi and F2 are 
equivalent to that of an effective ferromagnet F. Introduce the spin polarization vectors qi 
and q2, where g*^ has the magnitude of F^P^ and the direction of the magnetization direction 
of F^, with /3 = 1, 2. It is easy to test that these vectors obey the vector composition rule, i.e.. 
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q—qi+q2, in which ^is the spin polarization vector of F. Therefore, the effective parameters 
of F are 

f = ri + r2 , (23) 

p ^ [{riPiY + (r2P2)^ + 2riPir2P2cosg]^ 
ri + r2 

As a result, the total conductance of Fi and F2 can be obtained as 

G = G, + G2^GFDs{P,f) , (24) 

in which Gpus has the same form as in Eq.(22), P is the effective polarization, and f 
is defined by V g/V . Then the conductance of Fi and F2 can be expressed by the total 
conductance multiplied by a sharing factor, 

^ _ ^ + rir2 - {vlPf + v,p,V2P2 cqs e) 

^ rf + ri + 2rir2-(r?p2 + r2p| + 2riPir2P2cos^) ' ^ > 

^ _ ^ + rir2 - (r^p| + riPir2P2 cos e) 

' T\ + Tl + 2T^T2 - {TlPf + riP| + 2riPir2P2 COS 9) - ^ ' 

Fig. 3 shows the curves of G vs ^ (also can be viewed as 2Gi vs 9 or 2G2 vs 9 ) for the 
symmetric case, in which Fi = r2 = F and Pi — P2 = P. For r — 1, G increases with the 
increase of 9 or decrease of P. For r > 1, the curves of G vs ^ is quahtatively the same 
as those of r = 1, but the conductance is lowered and more sensitive to P. For r < 1, the 
variation is more complicated: if P^ < 1 — r^, G decreases with the increase of 9 or decrease 
of P; if P^ > 1 — r^, G has the maximum ^ at 9 satisfying (Pcos|)^ = 1 — r^. These 
results are readily understood by the new matching condition P^ + = 1 with the effective 
spin polarization P = Pcos |. 

Two points are noteworthy in the above result: (1) If Fi and F2 are regarded as a 
whole, the effective polarization can be tuned continuously by changing the angle of the 
mutual orientations, which is impossible for one chosen ferromagnet. (2) For r ^ 1, the 
total conductance for the two ferromagnets ins anti-ferromagnetic alignment is larger than 
that in ferromagnetic alignment, which is completely different from the effect of GMR or 
TMR. To describe this new effect of magneto-resistance, define the ratio of Andreev reflected 
magnetic resistance (ARMR) in (Fi,F2)-QD-S by 

ARMR ^ . (27) 

and the curves of ARMR vs P for various r are shown in Fig.4. 
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Fig. 5 shows the curves of Gi vs 9 for an asymmetric case, in which Pi = 1, P2 is arbitrary, 
and Fi = r2 = ^s/"^- Since Fi is fully polarized, the conductance of Fi is sensitive to the 
spin polarization and orientation of F2. For P2 = 0, Gi does not depend on 9; while for 
P2 = 1, Gi strongly depends on 9, with Gi = at ^ = and Gi — ^ at 9 = n. We suggest 
that this effect can be applied to measure the spin polarization of F2. In practice,one may 
chose a half-metal material as Fi, the ferromagnetic material to be measured as F2, and 
changing the spin orientation of Fi by applying an external magnetic field, then the spin 
polarization of F2 can be deduced from the weak / strong dependence of Gi on 9. 



IV. FINITE BIAS CURRENT 



Now we turn to investigate the non-equihbrium transport of (Fi,F2)-QD-S. For simphcity, 
we only consider the antiparallel orientation of Fi and F2 (i.e., 9 = 71), with finite but small 

bias voltages (i.e., |eVi| < A and |eV2| < A). Notice that the self-energy becomes to 
block-diagonal due to 6* = tt, and the expression of current Ii can be simplified as. 



h = ^j du [An(/i - /i) + Ai2(/i - h) + Qisih - fs) + Qi2{h - /2)] , 



(28) 



^11 


— Ti^Fi-f 


1 z^** 
l*-^12 




1 r^l' 1 


A12 


— r2Tri-i- 


\^12 


' + r2iFi^ 


|/-ir 1 



Qls — ^ll^sP 



Q 



12 



r2^Fi| I 



\G[,\' +\G\,\^ + 2Re{--G\,Gr2) 



\Gl,\' +\Gl,\^ + 2Re{+-Gl,Gll) 



G\i\ + r2trii \G 



2TJ- U 



•^331 



At zero temperature and in the low bias regime, the current of Qu process vanishes. 
Further assuming that both Fi and F2 are fully polarized, both Q12 and An process are 
also forbidden. Only the process of A12, i.e., crossed AR involving Fi,F2 and S contributes 
to the current, /i and I2 are derived as 



(ia;F2|Fi| \G'[^f' (/i - /2) 



(29) 



Notice that /i = I2 holds even if Fi ^ F2 and Vi ^ V2, because /i is pure spin f current 
and I2 is pure spin [ current while /i -|- I2 is required to be non-spin-polarized current by 
the superconductor. For simplicity, we further assume that Fi| = F21 = Fl (Fij^ = F2j^ = 
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due to Pi = P2 = 1) and Ts = Tr, then the system (Fi,F2)-QD-S is similar to a special 
N-QD-S one, in which the two spin bands of N have different chemical potentials controlled 
by Vi and V2. Define the transmission probability of crossed AR by Tar{uj) = ri|G^2r,the 
current can be expressed as 



e r^' 

I=T TAR{uj)du . (30) 
n .I_V2 



Notice that Tar{uj) is an even function of u, the above current formula implies that the 
sign of Ji or I2 is not determined by Vi or V2 but by |(V^i + V2). This is quite unusual because 
it contains the case that Vi > and V2 < but Ji = /2 > (This unusual property was first 



addressed in pO[ ) . Generally, for a three-terminal system, one may expect that current flows 
out of the terminal with highest voltage and into the terminal with lowest voltage. But for 
the current conducted by crossed AR, the sign of current in each ferromagnetic terminal is 
linked to the averaged chemical potential of the two, because two ferromagnets cooperate 
with each other in this process, with total energy balanced. Fig. 6 illustrates the conducting 
process corresponding to the case of Ji = /2 > with /ii > and /i2 < but |(/ii + 1x2) > 0. 

We ignore the energy structure of QD in Fig. 6 for simplicity, however, the current 
I = Ii = I2 depends strongly on the transmission probability of QD. In fact, I is the 
integral of Tar{uj) over the range of (— V2, Vi). Fig7 shows the surfaces of /(Vi, V2) and cor- 
responding Tar{uj) spectrum for three typical cases of Tl and Tr. In Fig. 7a, <^ Tr, the 
spin degenerate level of QD is hybridized to two Andreev bound states due to coupling with 
S, while the coupling with Fi and F2 provides the small broadening to these bound states. 
Tar has two peaks with the maximum of unity at each of the Andreev bound states. Corre- 
spondingly, the surface of /(Vi, V2) has five steps: the highest step corresponds to (— V2, Vi) 
covering both of the peaks; the second step (including two patches) corresponds to (— V2, Vi) 
covering one of the peaks; the third step (including three patches) corresponds to (— V2, Vi) 
or (Vi, —V2) covering none of the peaks; the fourth step (including two patches) corresponds 
to (Vi, — V2) covering one of the peaks; and the lowest step corresponds to (Fi, — V2) covering 
both of the peaks. In Fig. 7b, F^ = Vr, the Andreev bound states are sufficiently broadened 
so that the two peaks in Tar merge into one. The one peak structure of Tar spectrum cor- 
responds to three step pattern in /(Vi, V2) surface. In Fig. 7c, F^ F/j, the resonant level 
of QD is significantly broadened, as a result. Tar is small and flat with tails at u; = ±A. 
The structureless Tar spectrum corresponds to a plain in /(Vi, V2) surface, proportional to 
|(Vi + V2). In short. Tar spectrum can be extracted from the measurement of /(Vi, V2) 
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surface. 



V. CONCLUSIONS 



In this paper, we have investigated the Andreev reflection in a (Fi,F2)-QD-S system. 
By using the non-equihbrium Green function, a general current formula is derived, allow- 
ing arbitrary spin polarizations, magnetization orientations and bias voltages in Fi and F2. 
The formula is applied to several special cases, revealing some interesting properties of this 
system: (1) Analogous to the Fermi velocity mismatch in F/S contact, the zero bias con- 
ductance in F-QD-S reaches its maximum ^ if matching condition r^-fFij^ = F^j satisfied. 
(2) For total current (conductance) of (Fi,F2)-QD-S with Vi — V2, the two ferromagnets Fi 
and F2 are equivalent to an effective ferromagnet F, and the effective polarization P can be 
tuned by the angle between the spin orientations of Fi and F2. (3) There is a new effect of 
magneto-resistance in (Fi,F2)-QD-S (named as ARMR), in which the conductance for Fi 
and F2 in anti- ferromagnetic alignment is larger than that in ferromagnetic alignment. Base 
on this effect, a possible way to measure the spin polarization of ferromagnetic material is 
proposed. (4) The non-equilibrium transport of this system is quite unusual. Especially, if Fi 
and F2 are fully but opposite polarized, the signs of current through Fi and F2 is determined 
by ^{Vi + V2) rather than Vi or V2. Furthermore, the surface of I{Vi, V2) depends strongly 
on the AR transmission probability, which can be applied to extract the latter. Finally, we 
believe that the suggested (Fi,F2)-QD-S system is accessible of the up-date nano-technology, 
and we are eager to see relevant experiment on such appealing system. 
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FIGURE CAPTIONS 



Fig. 1 Schematic drawing of the three-terminal system under consideration. Fi and F2 rep- 
resent two ferromagnetic electrodes with different magnetization orientations and bias 
voltages, QD is a quantum dot, and S is a superconductor with zero voltage as the 
ground. 

Fig. 2 The zero bias conductance G vs P for F-QD-S, where P is the spin polarization of F. 
r = Tji/Tl is the ratio of coupling strengths, with r ^ 1 for (a) and r ^ 1 for (b). 

Fig. 3 The total conductance G vs 6 for (Fi,F2)-QD-S, where 9 is the angle between the 
orientations of Fi and F2, with Fi = r2 = F, Pi — P2 = P, and r = ro/(ri -|- r2). 

Fig. 4 ARMR vs P in (Fi,F2)-QD-S, where ARMR = {Gaf - Gf)/{Gaf + Gp), P and r 
have the same meaning as in Fig. 3. 

Fig. 5 The conductance Gi vs 9 for different P2, with P^ — 1 and r — 1. Gi has strong / 
weak dependence on 9 for large / small P2, which can be applied to measure the spin 
polarization of F2. 

Fig. 6 Schematic diagram of non-equihbrium transport in (Fi,F2)-QD-S. Fi and F2 are in 
anti-ferromagnetic alignment, marked by left- and right- slanted shadows, respectively. 
S is marked by crossed shadow, with the energy gap region ±A with respect to the 
chemical potential; QD is between the two barriers, and the energy structure is ignored 
for simplicity. The diagram illustrates an unusual property of the current conducted 
by crossed AR involving Fi, F2 and S: the signs of /i or I2 are determined by ^{1^1+^2) 
rather than /Xi or 112- 

Fig. 7 The current I = Ii = I2 vs the bias voltages (Vi, V2) for three typical cases: (a) ^ 
r^; (b) Tl = Tfi] and (c) ^ Fr. The surface of /(Vi, V2) has close relationship to 
the spectrum Tar{u!), which can be used to extract the latter. 
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Fig.4 



(a) rL=0.02A rp=0.5A 




(b) rL=0.lA rp=o.lA 





Fig.7 



